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Abstract 

Gene delivery of nucleic acid to the cell nucleus is a fundamental step in gene 
therapy. In this review of modeling drug and gene delivery, we focus on the particular 
stage of plasmid DNA or virus cytoplasmic trafficking. A challenging problem is to 
quantify the success of this limiting stage. We present some models and simulations 
of plasmid trafficking and of the limiting phase of DNA-polycation escape from an 
endosome and discuss virus cytoplasmic trafficking. The models can be used to assess 
the success of viral escape from endosomes, to quantify the early step of viral-cell 
infection, and to propose new simulation tools for designing new hybrid-viruses as 
synthetic vectors. 

Keywords: Mathematical Modeling, Gene Delivery, Stochastic Processes, Cytoplasmic 
Trafficking, modeling early steps of infection, Endosomal Escape, Narrow escape. 

1 Introduction 

Drug delivery is a multi-step process of delivering a cocktail of molecules through different 
tissue levels to a specific location in the body, which is usually a difficult task. For example, 
the Brain Blood Barrier is a physical barrier that drastically limits the access of molecules 
to the central nervous system, preventing efficient drug delivery. The possibility of injecting 
to a cell a piece of DNA that can reach the nucleus, in order to synthesize a given protein, 
is a very seducing idea, which has driven gene therapy and genetic crop modification. Gene 
therapy has been tried in the quest to cure diseases caused by single-gene defects, such as 
cancer and hereditary diseases linked to a genetic defect [1]. 

Non-viral methods for gene delivery use a variety of tools such as polymeric gene carriers, 
microinjection, gene gun, hydrostatic pressure, electroporation, continuous infusion, and 
many others [1]. Gene delivery can be used, in principle, for fighting major diseases, ranging 
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from viral infection to epilepsy. Although viruses are great predators, they can serve as 
vectors of drugs. For example, the Adeno-Associated Virus of type 2, (parvoviruses 
can carry gene coding for the neuropeptide Y, which is transduced in some areas of the 
hippocampus, where it is injected. The peptide is then released from neurons during high 
brain stimulations and reduces glutamate release, thus preventing the propagation of a high 
firing rate in neighboring neurons This therapy has the potential to reduce epileptic 
form activity. 

In this review, we discuss recent models of intracellular delivery and viral trafficking. 
Mathematical and physical models are constructed for the purpose of predicting and quanti- 
fying infectivity and the success of gene delivery. The models give rise to rational Brownian 
dynamics simulations for the study of sensitivity to parameters and, eventually, for testing 
the increase or the drop in infectivity, by using simultaneously a combination of various 
drugs. The modeling approach can be used for the optimization of the delivery in a high- 
dimensional parameter space. In the first part of the review, we present several approaches 
to the study of some limiting steps of gene delivery. A major limiting step remains the 
success of plasmid escape from an endosomal compartment. In the second part, we present 
models of viral entry, the endosomal step and trafficking in the cytoplasm. These models 
can be used to predict the efficacy of new synthetic vectors, based on hybrid viruses. 

2 Mechanisms of gene delivery. 

Introduction: Two types of gene vectors 

Gene therapy is a new therapeutic strategy that offers the promise of treating diseases 
through insertion, alteration, or removal of genes within an individual's cells and biological 
tissues. In general, any drug molecule must reach its intended site of action to exert its 
effect and to avoid non specific potentially dangerous interactions. A daunting hurdle of 
drug delivery is the large amount of charged molecules such as DNA that need to reach 
the cell nucleus. Indeed, the DNA molecule must initially pass through the cell membrane, 
traffic inside the cell cytoplasm to finally enter the nuclear pore in a form suitable for 
transcription to start. To overcome efficiently these barriers, the genetic material is usually 
associated to a viral or synthetic vector (see figured]). Viruses have developed evolutionary 
tools to enter cells and get transported towards the different compartments such as the 
nucleus and finally to reproduce using the protein synthesis machinery of the host cell (see 
figure dJ-b). However their use presents several limitations: viral vectors are not safe and 
have been implicated in death during clinical trials |4j, they can trigger an immune response 
and thus cannot be administrated repeatedly. Finally, the size of the transferred genome is 
limited by the nucleocapsid architecture. To overcome the difficulties imposed by the use 
of viral particles as vectors, synthetic vectors such as lipids, cationic polymers, peptides or 
combinations thereof, have been developed in parallel. In particular, the most common used 
vectors for nonviral gene delivery are cationic lipids, first introduced in 1987 [5] and later 
on cationic polymers [6j, followed by improved [7j polycations such as the polyethylenimine 
(PEI) for transfection experiments. Cationic lipids are macromolecules consisting of an 
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vector (PEI) 



Figure 1: Viral and synthetic vectors are used to transfer genes in the nucleus, (a) 

characteristics of viral and synthetic vectors. While synthetic vector polycations are preferred to 
viruses, they are safer and the size of the genes is not limited by any capsid architecture, but 
their efficiency to escape endosomes and traffic through the risky cytoplasm is still very low. (b) 
schematic description of early step of infection for viral and synthetic vectors. Synthetic vectors 
are not assisted by active transport during their cytoplasmic trafficking. 

aqueous core enclosed in a spherical phospholipid layer with positively charged head groups, 
while polycations are positively charged polymers which can have different structures (linear 
or branched). Both molecules are able to condense and neutralize DNA molecules through 
electrostatic interactions between their positive charges and the negatively charged DNA 
phosphate groups. In both cases, polycations protect DNA from undesirable degradation 
and facilitate entry into cells. 

Despite the high degree of reproducibility, the ease of modification, the lack of toxicity 
and the ability to deliver large pieces of DNA, polycations and liposomes are still much 
less efficient than their viral counterpart due to low cellular uptake and endosomal escape. 
A new promising strategy is the use of cationic liposomes incorporating fusogenic peptides 
from viral glycoproteins, which are supposed to enhance endosomal escape by mimicking the 
mechanism of fusion of viral envelopes with host cell endosomal membrane ^ |^. 

Unraveling the molecular mechanisms that underly the cellular behavior of both viruses 
and synthetic vectors is now needed to develop and optimize efficient hybrid vectors. In 
particular, the endosomal escape remains a major barrier in gene delivery and prompt re- 
lease from an endosomal compartment presumably constitutes one of the critical steps in 
determining the efficiency of transfection. In addition, the cytosolic motion of large DNA 
molecules is limited by physical and chemical barriers of the crowded cytoplasm [ini ITT] : 
whereas molecules smaller than 500 kDa can diffuse, larger cargos such as viruses or syn- 
thetic gene vectors, require an active transport system [12] such as the microtubules (MTs). 
Consequently, to understand how molecular components of gene vectors can affect quantita- 
tively their ability to pass through these two main limiting barriers of gene expression, we 
will present a review on our recent biophysical models for both viral and synthetic vectors. 
Indeed, while modeling the endosomal escape of viruses will help to understand the molecular 



3 



mechanisms underlying their rehable escape at a given pH, a model of the free cytoplasmic 
step will allow to analyze the active transport along MTs, which increases the probability 
that a viral particle reaches a nuclear pore in comparison with a Brownian synthetic gene 
vector. 

The cellular uptake pathways 

Since the phospholipid bilayer of the plasma membrane has a hydrophobic exterior and a 
hydrophobic interior, any polar molecules, including DNAs and proteins, are unable to freely 
pass through cell membranes. To circumvent this barrier, an approach consists in transfer- 
ring naked nucleic acids by inducing a membrane destabilization. An externally applied 
electrical field can temporarily disrupt areas of the membrane allowing polar molecules to 
pass, then the membrane can reseal quickly and leave the cell intact fi3\ [HI [15]. This 
method, commonly known as electropermeabilization or electroporation, provides versatility 
and efficiency, but it suffers from non specific transport across the membrane which can lead 
to ion imbalance and thus improper cell functions. It can also lead to irreversible cell damage 
if the electrical field is too strong. Microinjection also permits rapid delivery of genes to 
the cytosol or the nucleus by using a glass micropipette to inject genetic material into cells. 

Permeabilization is another used technique of delivery in which pore-forming agents which 
have the ability to fuse with the membrane are used to form large apertures in the cell 
membrane. However these techniques are highly invasive and cannot be used for in vivo gene 
delivery. A different strategy to overcome the initial cell membrane barrier is by associating 
genetic molecules with viral or synthetic vectors such as hpids, cationic polymers, peptides 
or combinations thereof, which are internalized through the endocytic pathway. Endocytosis 
is the process by which cells absorb molecules by engulfing and enclosing them into vesicles. 
There are several different endocytic pathways: 

1. Clathrin- Mediated Endocytosis Clathrin-Mediated Endocytosis (CME) is a process by 
which cells internalize molecules through the strong binding of a ligand to a specific 
cell surface receptor. This process results in the clustering of the ligand-complexes in 
coated pits which are formed by cytosolic proteins, the main unit being clathrin. The 
coated pits then invaginate the plasma membrane to form clathrin-coated vesicles. 

2. Caveolae-Mediated Endocytosis Caveolae are small (approx. 50 nm in diameter) pits 
in the membrane that resemble the shape of a cave which can mediate uptake of extra- 
cellular molecules through complex signaling. 

3. Macropinocytosis A mechanism of endocytosis in which large droplets of fiuid are 
trapped underneath extensions (ruffles) of the cell surface. 

4. Phagocytosis is the process by which cells bind and internalize large (> 0.5yu) pathogens 
such as bacteria or micro-organisms. 
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Endosomal escape 



Once a gene vector (synthetic or viral) enters an endosome, it has to escape into the cyto- 
plasm before being degraded in lysosomes. Although the exact pathways leading to endo- 
somal escape are not fully elucidated, they are hmiting steps in gene dehvery. Most viruses 
possess efficient endosomolytic proteins allowing them to disrupt the endosomal membrane, 
such as the VPl penetration protein of the adeno-associated virus (AAV) [T6] or the in- 
fluenza hemagglutinin (HA), [T7]. In addition, the biophysical mechanism leading to endo- 
somal membrane destabilization and concomitant plasmids release for synthetic vectors is 
still poorly understood. However, in both cases, acidification of the endosome is needed to 
trigger endosomal escape. For viruses protons or low pH activated proteases bind viral en- 
dosomolytic proteins, triggering their conformational change into a fusogenic state pTf [T6] . 
Polyplexes (DNA/cationic polymer) use a different strategy to escape endosomal compart- 
ments, based on the so-called proton sponge mechanism. Indeed some polycations have a 
high buffering capacity in a wide range of pH thus preventing the DNA from degradation 
and inducing an influx of chloride CI" ions to maintain electrical neutrality. The osmotic 
pressure due to the influx oi Cl~ ions causes endosome swelling and eventually membrane 
disruption. 

Cytoplasmic trafficking 

Following the endosomal escape, viral and synthetic gene vectors have to travel through the 
crowded cytoplasm to reach the nucleus and deliver their genetic material through the nuclear 
pores. While the cytoplasmic movement of viral particles towards the nucleus is facilitated by 
the microtubular network and viral proteins, very little is known about the fate of non-viral 
DNA vectors in the cytoplasm. However, trapping of large DNA particles (>500 kDa) in the 
crowded cytoplasm drastically hinders their cytoplasmic diffusion [101 E] and subsequently 
diminishes the transfection rate of synthetic gene vectors. In addition, one of the critical 
barriers in polyplex-mediated gene delivery is the timely unpacking process of the complexes 
within the target cell to liberate the DNA for efficient gene transfer [18]. A desirable property 
of a gene carrier is that it strongly binds to DNA to form more compact structures that can 
provide better protection of genetic material against nucleases and efficient transportation 
through the cytoplasm. However stronger polycation-DNA interactions produce a counter- 
productive effect on the timely release of DNA for transcription. Therefore an ideal vector 
would be capable of binding strongly to the DNA during the early stages of intracellular 
transport and to release the DNA right before nuclear entry. 

Nuclear delivery 

Expression of therapeutic proteins in gene therapy requires that the transfection agent, a 
virus or a synthetic vector, delivers its cargo through the nuclear envelope (NE) into the 
cell nucleus. Later on, the genetic material will be involved in the cellular transcriptional 
machinery. Some viruses such as AAV are able to transiently disrupt the NE delivering 
their intact capsid for disassembly inside the nucleus. However, most viruses translocate 
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inside the nucleus through nuclear pore complexes (NPCs) by using nuclear localization se- 
quences (NLSs). This facilitated nuclear import can accommodate the transport of molecules 
with diameters of up to 39 nm [2D112T]. Thus, small viruses such as hepatitis B virus [22] or 
baculoviruses |23] can cross the NPC and disassemble inside the nucleus. In contrast, larger 
viruses such as HIV [21] or influenza [25] uncoat inside the cytoplasm, or near a NPC, for the 
adenovirus p6j. The released components contain NLSs and are thereby able to cross the 
NPC. Viral derived NLS peptides associated to polycations have also appeared to enhance 
the nuclear delivery of genes |27j . 

Recent biophysical models have been proposed to quantify the molecular mechanisms 
of viral facilitated transport through the NPC |28l |29l |30l [31], however the exact import 
mechanisms still remain unclear. Interestingly, it has been found that pressure due to DNA 
condensation in double-stranded DNA viruses, such as the herpes virus, is as high as 50 
atmospheres, allowing a direct injection of the genome inside the nucleus through the NPC 



3 Modeling the endosomal pathway in Gene Dehvery 
with cationic lipids and polycations. 

The relative inefficiency of transfection using polycations, compared to the use of viral vectors 
remains the largest barrier to synthetic vector development and applications. Elucidating 
how non-viral vectors behave at the intracellular level is therefore a necessary step for syn- 
thetic vector improvement and optimization [33]. Biophysical modeling of polyelectolytes in 
confined domains, their interaction with biological membranes and escape from endosomes 
would allow us to design new efficient vectors. Computational models have accounted for 
endosomal escape using kinetic equations [31], where parameters describing multiples steps 
of gene delivery pathways are specified from experiments, numerical fits or through statisti- 
cal mechanics analysis [35]. However, it is difficult to extend these methods to identify the 
biophysical mechanisms underlying DNA endosomal escape. 



Two biophysical scenarios to study the endosomal release of DNA 

Cationic polymers, such as PEI, condense DNA into nano-sized polymer/DNA complexes 
(polyplexes) by a self-assembling process due to electrostatic interactions of the positively 
charged polymer with the negatively charged DNA. Polyplexes with positive surface charges 
are formed when the number of positive charges of the polymer exceed that of the neg- 
ative DNA charges. When a polyplex encounters the cell surface, it interacts with the 
negatively-charged cellular membrane and can be taken up into the cells via endocytosis. In 
the intracellular environment the polyplexes are located in endosomes that become acidified. 
In this case, DNA is prone to degradation by lysosomal enzymes. In order to transfer DNA 
cargo to the nucleus, polyplexes must escape from endosomes. This transmembrane mecha- 
nism or endosomolytic process remains unclear and it is believed that high gene transfection 
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Figure 2: The proton sponge hypothesis: protons and chloride ions Cl~ enter into the 
endosome. The increased concentration in the endosome with respect to the cytosol causes an 
osmotic pressure on the endosomal membrane which can induce its disruption. Swehing mechanism: 
Protons bind to the polycation thus increasing its charge density and stiffness due to higher Coulomb 
repulsion between monomers. In the figure, the darker chain represents the DNA while the lighter 
one represents the polycation molecule. 

efficiency observed with PEI is due to its ability to avoid endosome acidification and to pro- 
mote complexe's release through a mechanism called "proton sponge". This mechanism can 
be summarized as follow ([2]): the PEI contains amine groups, which can be protonated at 
low pH. Thus due to an infiux of chloride ions Cl~ to maintain the overall electroneutrality, 
and water molecules, the endosome starts swelling leading to an increased osmotic pressure, 
a destabilization of the endosome and ultimately disruption. It has not been proved yet 
whether the stress produced by the proton sponge effect and by the swelling are the only 
mechanisms responsible for the endosomal membrane disruption. However, recent exper- 
imental evidences argue in favor of this escape mechanism [36]. There are other possible 
mechanisms responsible for endosomal disruption such as swelling of the endosome due to 
internal charge repulsion triggered by an increased of ion concentration in the endosome. 
This effect implies a direct physical interaction between the DNA-PEI complex and the en- 
dosomal membrane. The PEI-DNA complex can escape an endosome through a membrane 
hole, probably due to direct interaction of polymers with membrane |37]. For example, the 
interaction between the cationic polymer PEI and cellular components [381 EH] lead to Gram- 
negative bacteria membrane disruption by free PEI, |10]. Furthermore, PEI also causes a 
liposomal membrane permeability to increase ^41j. Thus, PEI is capable of destabilizing lyso- 
somes and this suggests that it enables the DNA conjugated to the PEI to escape into the 
cytoplasm. Finally, a wide variety of nanoparticles, including different cationic polymers can 
physically disrupt lipid membranes by forming nanoscale hole and membrane thinning [12] . 
Cationic polymers contribute to expand existing defects or can directly induce small holes. 
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Today, most common theories of pore formation are derived from the classical nucleation the- 
ory [53] , which unfortunately does not predict stable pores of finite size. However long-lived 
pores that remained open for several seconds have been observed [lH SSI SSI SZl SB] • Several 
theoretical models have succeeded to show that after pores nucleation, their enlargement 
is expected to relax the surface tension, leading to the formation of a stable or long-lived 
metastable pore gOl [501 EH [521 [53l El] . 

For a planar membrane, increasing the lipid density occurs concurrently to pore dilation, 
that reduces mechanical tensions. After a vesicle opens, its internal content escapes, reducing 
the osmotic pressure and the associated surface (Laplace) tension. Although membrane fluc- 
tuation was recently modeled [551 ISHl |57] , there is still no consistent theory of pore formation 
and in particular, it has yet to be elucidated whether polymer-membrane interactions and 
hole formation are responsible for endosomal escape and how such process can be controlled 
to optimize gene cytoplasmic delivery. 

Molecular simulation of DNA complexation in an endosome. 

To check whether protonation increased of a PEI chain is associated with endosomal swelling 
disruption, we summarize here some recent findings [SB] about complexation between two 
oppositely charged polymers. Other theoretical studies [5S1 [SD] have shown that after com- 
plexation, extended oppositely charged polymer chains collapse into a compact globule, even 
though the initial chain configuration contained bound ions. Ions are progressively released 
by polycation-polyanion pairs attraction and the final polymer complex, containing two 
chains separated from the rest of the ions. 

In the presence of many interacting ions, the dynamics of a DNA molecule mixed with 
polycations is quite complex. But motivated by experimental data [61] suggesting that 
the number of polycations can affect the PEI-DNA interaction, we recently used Brownian 
simulations |58] to study the effect of various proton concentrations on this interaction. 
We modeled polymers as flexible bead-spring chains, where each bead represents a charged 
monomer, while ions are modeled as charged particles. The solvent is treated as a dielectric 
continuum with dielectric a constant e = 80. We started with two oppositely charged chains 
whose monomer positions are respectively 



where the charge electro neutrality is preserved P = M + N. The motion of a polymer 
chain in an overdamped medium is described by the Smoluchowski's limit of the Langevin 
equation, which reduces to a system of first order stochastic differential equations. The 
dynamics of bead i at position Xi in a potential U{x, y, z) is 




(1) 
(2) 




mixed with ions at position 




(3) 




(4) 
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where D is the diffusion coefficient and tbi, ..tiyjv are N- independent 3— dimensional Brownian 
motions and 7 is the viscosity coefficient. The total potential U — U{x, y, z) is the sum of 
various potential terms 

U (x, y, z) = Uel{x) + UEL{y) + Ubend(x) + Ubendiv) + Ulj(x, y, z) + Uc{x, y, z) (5) 

where 



• Uel{x) is the polymer elastic energy: the potential well Uk generated by two adjacent 
springs on the A;-bead is the sum of the two neighboring potentials: 

Ukk+i{xk, Xk+i) = k (^-^\xk - Xk+i\'^ - lo\xk - Xk+i\ 

Ukk-i{xk, Xk-i) = k (^\xk - Xk-i\'^ - lo\xk - Xk-i\^ , 

where k is the linear elasticity constant and Iq the equilibrium length. For < k < N, 
the potential Uk is 

Uk{xk-i,Xk,Xk+i) = Ukk+i{xk,Xk+i) + Ukk-i{xk,Xk-i) 

Uo{Xi,X2) = Ui2{Xi,X2) 

Un{xn-i,Xn) = Unn-i{xn,Xn-i) 
and the total potential is given by 

N 

Uel{x) = ^ Uk{xk-i, Xk, Xk+i) (6) 

k=l 

The elastic energy UEhiv) of the second chain is computed similarly. 

• Ulj{x) is the repulsive Lennard- Jones (LJ) potential, which models the excluded vol- 
ume between all the particles in the system (monomers and ions) and it is computed 
as follows: the LJ-potential between two monomers of the same chain is given by 



ULj{Xi, Xj) = 



4eLj 



+ ^LJ Xij < X 



(7) 



*^ c 



where the distance between the monomers, e is the depth of the 

potential well, a is the (finite) distance at which the inter-particle potential is zero. 
The LJ-potential is cut off at a distance Xc — 2^/^(7 and the constant e^j is added to 
avoid discontinuity at Xc- 

A similar potential is used to compute the LJ-interaction between monomers of the sec- 
ond chain ULj{yi,yj), monomers belonging to different chains ULj{xi,yj), monomers 
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and ions, Uij^Xi, zj) and Uijijj^, Zj), ion-ion Uu^Zi, Zj), so that the total LJ potential 
is given by: 

Ulj{x, y,z) = ^ ULj{Xi, Xj) + ^ ULj{Xi, Vj) + ^ Ulj{ 

• The stiffness of a polymer is described by a bending energy of the form 

N N-1 
Ubend{x) = - 2Xi + Xi+iY = Kang ^(1 - • Wj+l), (10) 

i=2 i=l 

where Ui = (a?j+i — Xi)/\xi^i — Xi\ is the unit bond vector connecting two consecutive 
monomers and Xi is the position of the i — th monomer. This potential depends on the 
angle 9i between two successive monomers since Ui ■ Uj+i = cos9i . The constant Kang 
measures the bending rigidity which is related to the persistence length of the polymer 
by Ip = Kang/ {ksT), wherc ks is the Boltzmann constant and T is the temperature 
[62j . The persistence length characterize the stiffness of a polymer and it is defined as 
follows: 

{cosO) ~ exp — (11) 
Ip 

where 6 is the angle between two vectors tangent to the polymer separated by a distance 
/, and the cosines of the angles is averaged (angled brackets) over all configurations. 
The stiffness energy Ubendiv) of the second chain is similarly computed. 

• The electrostatic potential between two charged monomers (charges are relocalized at 
the monomer) at position Xi and Xj is given by the Coulomb potential 

t/cN.x.)^5^E^ (12) 

where e is the dielectric constant of the medium and Zi the valence of the ion i. The 
total electrostatic potential is 

Uc{x, y,z) = ^ Uc{xi, Xj) + ^ Uc{xi, Vj) + ^ Uc{xi, Zj) + (13) 

which represents the total Coulomb interactions between monomers on a single chain, 
Uc{Xi,Xj) and Uc{yi,yj)- For monomers located on different chains the potential is 
Uc{xi, i/j), for monomer- ion interaction, it is UciUi, Zj) and UciVi, Zj), and finally the 
ion- ion interaction is given by Uc{zi, Zj). 
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The dynamics of the second polymer with M beads and positions y = {y^,--- ^y^) is 
governed by 

y, + -— = V2Dw, for i = 1..M. (15) 
^dyi 

Finally the ions with positions z = {zi, ■ ■ ■ , zm) satisfy 

Zi + -— = V2Dwi for i = 1.. p. (16) 

lOZ, 



The previous stochastic equations were solved numerically [5B], by standard Euler's scheme. 
The results of the simulations are presented in the next section. 



Numerical simulations of DNA-polymer interactions 

To check whether increasing the proton concentration modulates the DNA-PEI interaction 
and leads to the deformation of the DNA-polycation complex, we ran various Brownian 
simulations [58], that we summarize now. We initially placed two oppositely charged chains 
of total length Lp in a 3D-cube of length L and each monomer position is randomly chosen 
at equilibrium (the two spherical angles between monomers are uniformly distributed in the 
2-sphere 5(2) of radius Iq). To ensure electroneutrality, we added a number of ions, equal to 
the difference between the DNA and polycation charges. Periodic boundary conditions are 
imposed on the cube. L was chosen such that the density p = {Lp/L)^ ^ 10~^ is small. In 
that regime, there are no self-interactions between polymers. In a first transient regime, we 
waited for each chain to relax to equilibrium. Then, we move the chains close enough and 
restart the simulations to evaluate their interaction. We found that the structure of the final 
complex made of the two chains, after equilibrium is that of a flexible polycation winding 
around the DNA chain (FigHj) and ions are uniformly distributed over the simulation box. 
To further study the consequences of decreasing the pH and changing the PEI ionization 
(accounting for the experimental procedure [SS]), we repeated the previous simulations in 
the presence of different ionic concentrations (Fig. Sj). 

We further run additional simulations in a reflecting sphere, modeling an endosomal 
compartment to estimate the distribution of charges along the PEI monomers. We used 
[58] a Metropolis Monte-Carlo simulation [61] to evaluate the protonation configuration of 
a polymer having N monomers. Each monomer is capable of binding a proton. To estimate 
this distribution, we compute the free energy at equilibrium by finding the minimum of 

^ 1 
G(x) = {x^kBT\og 10{pH - pK^)) + - ^ W,,q,q„ (17) 

where x = {xi,X2,--- ,xn} and Xi = 1 when site i is protonated, otherwise. The first 
term in[T7]is the energy required to protonate a site, depending on e, the pH and the pK of 
the PEI ^64] and Wij is the electrostatic Coulomb interaction energy between sites i and j 
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when they are both charged. Thus, the average protonation of a given site Xj is obtained by 
averaging (Boltzmann weighted sum) over all configurations: 



We used [58] the Metropolis algorithm [65] to sample the most probable states of Xi and 
Figj3] summarizes the configuration distributions for various pH. At physiological concen- 
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Figure 3: Simulated protonation curves of PEI for different salt concentrations, 

0.01M,0.05M, l.OM. At high salt concentration, there is a steep rise in the protonation as the 
pH approaches pK^ and the chain is fully protonated for pH — pK^ ~ 2. As the salt concentration 
decreases, protonation of the chain becomes more difficult, because there are fewer ions to screen 
the electrostatic repulsion of charges on the chain. 



tration Cg = 0.15M and pH = 7.4, then pH — pK ~ 2.7, that is 50% of monomers are 
protonated, higher than experimental measurements (20% of protonated sites |;66J). DNA- 
PEI complex seems to evolve to different configurations depending on the pH. In Fig. Hti, the 
DNA-polycation forms a complex. The polycation is only partially charged (indeed 20% of 
the amine groups are protonated [66] )• The polycation behaves almost like a flexible chain, 
while DNA monomers are subject to strongest Coulomb repulsion making it stiffer and more 
extended, and almost immobile during complexation. Fig. |4]d and Hb show the DNA/PEI 
complex conflgurations respectively at pH = 6 and 5, during the endosomal acidiflcation. 
Increased protonation could induce swelling [63] of the complex and consequently promoting 
endosome rupture. The result in |58] predicts that increasing attraction between oppositely 
charged polyelectrolytes lead to a more compact complex structure. Upon increasing the pro- 
tonation state of the PEI, the positively charged chain condenses more tightly on the DNA, 
which become more flexible because its charged are now screened. The resulting complex 
resembles a typical toroid conflguration found for DNA interacting with multivalent ions [67] . 
The same polyplex structure was previously obtained [68] by solving numerically modifled 
Poisson-Boltzmann equations for the electrostatic potential and the polymer concentration. 
These results [58] suggest that two oppositely charged polymers condense upon lowering pH 
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Figure 4: Snapshots of the equiUbrium configurations of the DNA/PEI complex at pH 

= 7.0 (left), pH = 6.0 (middle), ph = 5.5 (right). During the acidification process the size of the 
complex decreases. Each image represents the configuration at fixed pH after 10^ — 10^ steps of the 
Brownian dynamics, which corresponds roughly to 1 — 10s. 

and thus this mechanism cannot be responsible for the sweUing of the endosome and endo- 
somal membrane disruption. We shall now review another simulation sets, aimed to better 
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Figure 5: Snapshots of the equilibrium configurations for a DNA molecule complex with many 
small polycations, pH = 7.0 (left), pH = 6.0 (middle), ph =5.5 (right). During the acidification 
process, the complex-DNA can swell due to an increased Coulomb repulsion and consequently some 
polycations can be released from the complex. Each frame represents the configuration at fixed pH 
after 10^ — 10'' steps of the Brownian dynamics. 

characterize an experimental procedure where the DNA molecule is mixed in solution with a 
fixed concentration of PEI, where the ratio r = N/ P of the N amine to P phosphate groups 
is fixed. A value of N/ P = 4—10 has been reported from experiments [61] and we simulate 
a DNA molecule with Ni = 100 monomers when it is complexed with 20 polyactions, each 
of length A^2 = 20, making the ratio r = N/P = 4. We use the same numerical method 
as described above with many polycations. Fig. [5^ shows the equilibrium configuration of 
the DNA complexed with the small polycations at physiological pH with periodic boundary 
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conditions, while FigJSjj and|5t represent the influence of acidiflcation in a conflned spherical 
micro domain. At physiological pH, the DNA complex is less compacted and more flexible 
due to an excess of positive charges resulting in additional screening, compared with the 
situation of two chains only (FigH^). During acidiflcation, at higher charge fraction of the 
PEI, some polycations are released from the complex and are free to diffuse in the endo- 
some. The presence of free polycations has been also observed experimentally [6^. Recent 
experiments [TD] have shown that purifled polyplexes without free PEIs were less efficient in 
transfection compared to non-purifled polyplexes. Hence interaction of free PEIs with the 
endosomal membrane, could play a role in endosomal escape, however a detailed model of 
membrane destabilizaton due to membrane-polymer interaction is currently missing. 

In summary, Brownian simulations can be used to study processes occurring at interme- 
diate scale between molecular and cellular. We have reviewed here some approaches method 
to simulate the interaction of oppositely charged complex polymers in endosomal following 
protonation. We reported that surprisingly, the protons and chloride endosomal influx which 
lead to an increase in the ionic concentrations and the enlargement of the polycation-DNA 
complex (due to internal charge repulsion), in the range of parameters such as pH G [5.5 — 7], 
the number of poycations from 1 to 20 [58] , seems not to have any direct role in the endosome 
disruption. This result suggests that more reflned scenario should be considered to study 
the direct PEI-DNA interaction during endosomal membrane disruption. 

4 Modeling the endosomal step of viral infection 

To better understand how viruses can be used as gene vectors, we now focus on one of the fun- 
damental step of early infection, where the viral particles travel inside an endosome. Indeed, 
most viruses enter cells in an endosomal compartment, after binding to speciflc membrane 
receptors. To undergo cytoplasmic or nuclear replication [HI [72], and to avoid degradation 
in acidic lysosomes, viruses must then successfully escape the endosome. Enveloped viruses, 
such as Influenza, contain membrane-associated glycoproteins mediate the fusion between 
the viral and endosomal membranes. In particular, acidiflcation of the endosome triggers the 
conformational change of the influenza hemagglutinins (HA) into a fusogenic state, leading 
to endosome-virus membranes fusion and genes release inside the cytoplasm. 

Developing a biophysical model of the influenza endosomal step offers a general framework 
to study the influence of several parameters such as the endosomal size, the number of 
viral particles or the molecular structure of glycoproteins on the efficiency of viral escape. 
In particular, fusogenic peptides derived from viral glycoproteins are increasingly used in 
cationic synthetic vectors [El [9], and the pH-sensitivity of fusogenic glycoproteins can now be 
tuned by modifying the electrostatic stability of the fusogenic complex [73] . These engineered 
glycoproteins shall serve to design efficient gene vectors and quantitative models will help 
optimizing glycoprotein molecular properties with respect to the endosomal escape efficacy 
of the vector. 

We shall review previous models [^ aimed to estimate the residence time of a viral 
particle inside an endosomal compartment. These models are based on considering the 
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accumulation of discrete proton binding events, leading to the conformational change of 
HAs. This is the limiting step of genes release in the cytoplasm and the model was built into 
two steps: first, by fixing the concentration of protons, we used a Markov jump analysis to 
estimate the mean time for protons to bind the HA binding sites until a threshold is reached, 
triggering its conformational change into a fusogenic state. Interestingly, our analysis allow 
to extract from the HA conformational change kinetics measured experimentally at different 
pH [73] the binding rates and the number of bound sites that are needed for the protein 
to change conformation. Second, coupling the pH-dependent conformational change of HA 
glycoproteins with a linear proton influx rate, we analyzed the endosomal escape dynamics 
of influenza viruses. We predicted [71] that the size of the endosome drastically impacts 
both the escape kinetics and pH, which reconciles different experimental observations: while 
a virus can escape from small endosomes (radius of 80nm) in the cell periphery at a pH ~ 6 
in about 10 minutes [75], it can also be routed towards the nuclear periphery, where escape 
from larger endosomes (radius of 400nm) is rapid (less than one minute) at pH 5 [77]. We 
shall now review in details the modeling of these two steps. 

Modeling the conformational change of glycoproteins 

To estimate the conformational change rate of a single HA glycoprotein at a given proton 
concentration c, we considered [73] that the protein changes conformation instantaneously 
when the number of bound sites reaches a critical threshold Ucru < ng, where Ug is the total 
number of HA binding sites. 

To follow the time dependent number of protonated sites, we use a stochastic analysis 
based on Markov jump processes [751 [ZSl EHl EI], where during time t and t + At, the 
amount of protonated sites X{t,c) can either increase with a probability r{X,c)At when 
a proton binds to a free site, decreases with probability l{X,c)At when a proton unbinds 
or remains unchanged with probability 1 — l{X,c)At — r{X,c)At. Using a scaled variable 
x{t, c) = eX{t, c) where e = l/ug and Ax = x(t + At, c) — x{t, c), the transition probabilities 
satisfy 

Pr{Ax = e\x(t, c) = x} = r{x, c)At, 
Pr{Ax = —e\x{t, c) = x} = l{x, c)At, 
Pr{A,T = 0|x(t, c) = x} = (1 — r{x, c) — l{x, c)) At. 

For a fixed proton concentration c, the transition probability function p{y,t\x,c) that the 
proportion x(t, c) of protonated sites is equal to y at time t, x{t, c) = y, given that a;(0, c) = x 
is solution of the backward master equation [7S1 EHl EDI IHT] : 

p{y,t\x,c) = p{y,t — At\x + e,c)r{x,c)At + p{y,t — At\x — e,c)l{x,c)At 

+ p{y,t-At\x,c){l-r{x,c)At-l{x,c)At), (19) 
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which can be expanded as (Kramers-Moyal expansion) 



dp ^ 

— = L^p = r{x,c)22^_i9xTp{y,t\x,c) 

n=l 

-e 



+ /(x,c)^^^(9.)Xy,t|x,c). (20) 



n=l 



The first time r a glycoprotein is filled up to a critical threshold Xcrit = ncru/ng is the first 
passage time for the bound protons x(t, c) to reach the level Xcrit- The mean first passage 
time t{x,c) is defined as the conditional expectation r(x, c) = E[T\x(t = 0, c) = x], and 
satisfies Ell [82]: 



L^t(x,c) = -1 for X in [0,Xcrit\, 

OT ( X Ci 

r(x, c) = for X = Xcrit and — — = for x = 0. 

CJJb 

We approximate the conformational change mean time with to{c) = t{xo{c),c), where < 
Xo(c) < Xcrit is the mean number of bound protons (concentration c), the leading order in 
e < 1 [THl [801 EH is 

ro{c)^C{e,c)\l- ^P^] I, (21) 



viXcriti c) 



where 



C(6,C) 



27r 



£[l/r){xo{c),c) 



r(a;o(c),c) 0(xcrit,c) 



and 



1 r 

u \i I \ 1 — / loq [lis, c) Iris, c)) ds 
, lix.c) rix.c) - I e ^ ^ ■' 'i ^ ■' " 

^/(x,c)/r(x,c) 

Formula ([2T1) relates the conformational change mean time of a single HA at a fixed protons 
concentration c with the binding/unbinding rates r(x, c) and /(x, c) 



r(x, c) = is:cn,(l - x), and /(x, c) = l{x) = Kusil - x)10"(^(^~^)+^\ (22) 

obtained from the mean number Xq{c) of HA protonated sites at different pHs ^83j. In 
addition, = 9 is the number of binding sites [83j, c the concentration of free protons in 
the endosome and K the proton binding rate to free binding sites. Using these transition 
rates in formula ( l2Tll . and comparing the theoretical conformational change mean time to(c) 
with experimental data [75], we found that K pa 7.5 * 10^L.mol~^s~^ and Xcrit ~ 0.7 [75]. 
Interestingly, the theoretical curve To(c) describes the entire range of experimental data 
set of the infiuenza hemagglutinin [7S], confirming the validity of our approach. Finally, 
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comparing the mean proton binding time r = 1/ (Kc) with the mean time for a proton to 
find a binding site by diffusion [SI], we found [Hj that: 



Td/r ^ 10 



-4 



(23) 



suggesting that the HA binding time is dominated by a high activation barrier, which guar- 
antees the stabihty of a conformational change, that cannot be easily triggered randomly 
(table 2 [75]). The theory predicts that the HA conformational changes occurs when roughly 
Xcritns ~ 6 HA binding sites are protonated. 

Modeling the endosomal escape of the influenza virus. 

During endosomal maturation, protons enter actively through V-ATPase pumps located in 
the endosomal membrane, leading to pH decrease and the conformational change of HAs 
into a fusogenic state [T7]. After fusion of the viral and endosomal membrane, the influenza 
genes can be released inside the cytoplasm. 

For a linear time dependent proton influx Xt, the the rate A is proportional to the number 
of pumps and thus to the endosomal surface [H]. Consequently, A may drastically increase as 
the endosome matures and increases its size by fusion of early endocytic vesicles into larger 
compartments [85]. As protons accumulate into the endosome, they can bind to influenza 
HAs, triggering their change of conformation when exactly 6 sites are protonated. 

The entire viral membrane in contact with the endosome seems to fuse before genes 
escape, as observed in electron microscopy images (flgure 5c [76]). Furthermore, another 
step of endosomal escape consists in the enlargement of the fusion pore, that should rely 
on the activation of additional HA located nearby the contact zone, between the virus and 
the endosome membranes [86]. In previous models [71], we accounted for this complex 
cooperative mechanism between activated HAs nearby the contact zone, by considering that 
genes are released in the cytoplasm when the total number of activated HAs among the viral 
envelope reaches a threshold < T < Nha, where Nha = 400 [ST] is the total number of 
HAs covering a single virus. Consequently, the escape release time is deflned by 



where HAo{t), HAi{t) . . . HAQ{t) are the number of HAs that have bound 0, 1 ... 6 protons 
at time t. The acidiflcation time course of an endosome containing an influenza virus is 
related to the number of free protons located in the endosome by P{t) = jVVoc{t) , where 
c(t) is the associated endosomal concentration at time t. Using the on-rate of a proton to 
a HA free binding site, r (x) = r {x,P{t)) /P{t) = Kns{l — x)/{J\fV) and the off-rate l{x) 
(!22|) . the kinetics equations are 



= mi{t\HAQ{t) > T}. 



(24) 



HAi + P 



HAo + P 



r(0/ns) 
l{l/n,) 
f{l/ns) 
l{2/ns) 



HA 



HA2 



HA5 + P 



f(5/ns) 



^HA^ 



(25) 
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and the associated mass action law leads to the following differential equation system 



- x-j:f(^)pmA,,{t) + j2i(i-)HA{t) 



dt 



j=0 ^ i=i 



dt \''^sj xn. 

"^^^^^'^ - ^r(^]HMt)~f(-]HMt)]p{t) + l(-]HA,{t)-l(-]HA{t) 



dHMt) _ rJ^^HMt)P(i), (26) 



dt \n, 

where the proton influx rate in the endosome is a linear function of time At. The initial 
conditions (t = 0) are a neutral medium pH = 7 and Pq = A/'VolO"^ ~ 1, HAo(t = 0) = 
Nha and HAi{t = 0) = for 1 < i < 6. 

The two critical unknown parameters of the model are the protons influx rate A and the 
mean number T of activated HAs, needed for the large fusion pore formation allowing gene 
release inside the cytoplasm. These parameters were obtained by solving numerically the 
system of equations [221 aiid by comparing the time t^, (formula f2M with the experimental 
mean escape time (~ 10 minutes) obtained for the virus-endosome fusion in Hela cells [76] : 
we found [71] 

T ^ 50%Nha = 200, and A ^ 3s~\ (27) 



Using modeling approaches to study various endosomal pathways 

Viral fusion and genes release from small endosomes (radius ~ 80 nm) have been observed in 
cell periphery |76| , starting only 10 minutes after viral entry. Although the endosomal pH in 
these small endocytic vesicles is not precisely known, it is believed to be ~ 6. However, most 
viruses are routed towards the nuclear periphery into larger endosomes (radius ~ 450nm), 
where the escape time decreases to 1 — 2 minutes and the associated escape pH to ~ 5 [77] 
(FIG. 111). 

To model the escape process in these two endosomal pathways, we developed a model [71] 
that accounts for the endosomal size, where the proton influx rate is scaled to the endosomal 
surface A(r) = (r/ro)^ A = (r/(80nm))^ 3s~^ [71]. Taking this scaling into account, we found 
by numerical resolution of equations (1261) that the escape pH and mean escape time are 
drastically decreased as the endosome becomes larger. When the threshold to fusion is fixed 
to T = 200 activated HAs [71], the escape pH is 5.8, and the escape time is 10 minutes for an 
endosomal radius of 80nm, while the escape time and the pH drop to 40s and 5 respectively 
for a radius of 450r;,m (Fig. [Hl^-c). 

In the general context of cytoplasmic trafficking, we proposed ^74j that, when viruses use 
the periphery pathway where they leave the endosome far away from the nucleus, at a pH 
6, they can further escape the rapid inactivation process occurring at pH < 5.4 |89j. In that 
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case, they still have to travel through the risky cytoplasm to a small nuclear pore to deliver 
their genetic material. However, for viruses using the other pathway leading directly to the 
perinuclear area, the escape endosomal pH drops to < 5, thus this exposure leads to possible 
HA degradation resulting in greatly diminishing the escape capacity. Consequently, it is not 
clear which pathway is more efficient for viral gene delivery. 




T. (Beconda) 



Figure 6: Schematic representation of the influenza virus endosomal pathways. Genes 
are released in the cytoplasm and the mean escape time and pH depends on the 
endosomal radius (a) After cell entry, the coated vesicle (CV) containing the virus can either be 
routed towards small static early endosome (SEE) in cell periphery [881 176j. where fusion occurs 
after 10 minutes at pH 6 [76] (orange pathway), or it can fuse with dynamic early endosomes (DEE) 
traveling along microtubules, to reach the perinuclear region and rapidly fuse from larger endosomes 
(green pathway) [771 EH]- In the second case, fusion time is decreased to 1 — 2 minutes and the 
escape pH drops to 5. (b) Time evolution of the number of activated HAs HAeit) (grey line) and 
the endosomal pH (solid line) for a small SEE with r = O.OSfim [76J (proton influx rate A = 3s~^). 
The pH is maintained above 5.8 for HAQ{t) < 50%n//^ = 200 (highlighted with blue dotted lines) 
leading to genes release in the cytoplasm, (c) For a larger perinuclear DEE with r = 0.45/um [85], 
the proton influx rate has been scaled to the endosome surface A = (45/8)^ 3s~^ = 95s~^, which 
leads to a rapid gene release in about 40 seconds, d-e) Engineered Stabilized HAs with an increases 
low-pH stability shall prevent viral vectors escape from peripherical endosomes at pH 6. 



Perspectives in developing simulations to optimize synthetic virus-like vectors 

For many enveloped viruses, conformational change of glycoproteins into a fusogenic confor- 
mation is triggered either by cumulative discrete binding of protons (e.g. class I, II and HI 
fusogenic glycoproteins [901 El]), or by low-pH activated proteases (e.g. Ebola GP protein 
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and the SARS coronavirus spike protein S [92l|93]). In addition, fusogenic glycoproteins are 
increasingly used in many synthetic drug delivery systems such as siRNA delivery systems 
[HI] or cationic synthetic gene vectors Thus, the generic model developed in [71] for the 

influenza virus can be extended to this large class of viruses and synthetic vectors. Further- 
more, the pH stability of influenza glycoproteins has been recently tuned by incorporating 
or removing electrostatic bonds into the fusogenic complex [73]. Consequently, developing 
biophysical models to quantify the endosomal escape kinetics and pH of viruses with respect 
to glycoproteins molecular properties is timely and could be further applied to optimize gene 
vector design. Indeed approaches developed in [71] can predict the endosomal time course 
of vectors or mutant viruses covered by these HA variants. In [7^, we have predicted that 
using stabilized forms of HA glycoproteins [73] shifts the escape pH in small endosomes from 
6 to 5, which may prevent vectors to use the peripherical endosomal pathway (flgure[6]d-e). 
Consequently, if genes degradation inside the cytoplasm is a limiting barrier to transfection 
efficiency, we predict that stabilized mutant viruses may increase the ability of vectors to 
deliver genes to nuclear pores. 

5 Modeling gene carrier cytoplasmic dynamics 

We end this review by a section on quantifying the success of cytoplasmic trafficking. Indeed, 
following the endosomal escape, viral and synthetic gene vectors have to travel through the 
crowded and risky cytoplasm to reach the nucleus and deliver their genome through the 
nuclear pores. Recent imaging techniques now allow single particle tracking [951 ED EE], 
which can be used to model cytoplasmic motion of gene vectors. Because gene vectors do 
not possess means of locomotion, they entirely rely on diffusion and active transport along 
MTs to reach the nucleus. 

The cell cytoplasm is a highly crowded environment and diffusion of macromolecules de- 
pends on their size: while non interacting spherical particles with radius up to ~ 2'bnm are 
freely diffusible in the cell cytoplasm [^, increasing the size above Abnm reduces consid- 
erably the motion [11]. Interestingly large viral particles have developed nuclear targeting 
signals to be actively transported along microtubules, which resulted in drastically decrease 
their arrival time to the nucleus [12]. It remains a challenging question to analyze their as- 
sociated viral trajectories, which consist of a succession of free or confined diffusion and/or 
ballistic periods [13 [TH [HS]- The analysis of such random trajectories starts with the posi- 
tion X(t) at time t of the gene vector, which is a stochastic process [8T1[98] and the dynamics 
depends on forces applied on the particle. When the motion of the gene vector is purely 
diffusive such as for cationic synthetic vectors, the overdamped equation for the velocity 
is simply cD^/dt = \/2D dW /dt, where W is the standard Brownian motion. For a viral 
gene vector, the motion is usually assisted by active transport along MTs and the stochastic 
equation for the velocity becomes 



where b is a drift that accounts for ballistic periods along MTs. This continuous Langevin 
description can be used to generate computer simulations of trajectories |8T1 ES] in free and 
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confined environment [100] , and is the basis to derive asymptotic formulas for the probability 
and the mean first passage time of the vector to a nuclear pore [SI] . However, to reduce the 
complexity of large simulations, a first step consists in studying the dependency of b (X) as 
a function of the MTs organization and the viral dynamical properties (diffusion constant 
D, affinity with microtubules and net velocity along MTs) \101\ 1102] . 



Prom the stochastic description of vector trajectories to the probabihty and the 
mean arrival time to a small nuclear pore. 

The crowded cytoplasm is a risky environment for gene vectors that can be either trapped or 
degraded through the cellular defense machinery. Consequently, the cytoplasmic trafficking 
is rate limiting for genes expression, and to analyze quantitatively that step, we derived 
[1031 1104] asymptotic expressions for the probability P„ and the mean time r„ a single gene 
vector arrives to one of the n small nuclear pores. In particular, obtaining these expressions 
allows to explore the phase space of parameters and in particular links global quantitative 
outputs measuring the success of nuclear genes delivery with the cellular geometry, the MTs 
organization and the dynamical properties of gene vectors. To obtain these expression, we 
modeled the viral degradation or immobilization by a steady state degradation rate /c(x). 
We then introduce the survival probability density function (SPDF) p{x,t), which is the 
probability to find a live (not degraded) viral particle inside a cytoplasmic volume element 
X + dx ai time t [103] 

p{x, t)dx = Fr{X{t) e x + dx,T^ > t,T'' > t\ pi}, (29) 

where r" is the first time for a live virus to arrive to an absorbing nuclear pore, denoted 
dNa and the first time that it is degraded. The viral initial distribution is Pi. The SPDF 
p{x, t) satisfies the Fokker-Planck equation (FPE) [81j 

^P^^'^^ = DAp{x, t) - V ■ b{x)p{x, t) - k{x)p{x, t) for x eQ 

p{x,Q) = Pi{x) for a; G fi, (30) 

which describes the time evolution of the vector probability density function. The boundary 
conditions are 

p{x,t) = for X E dNa 
J{x,t)-nx = xedn-dNa, (31) 

where the first condition describes the particle absorption at the nuclear pores area dNa 
and the second one the reflection of the vector on the remaining boundary area of the cell 
dQ — dNa- The flux density vector J{x,t) is defined by 

J{x, t) = -DVp{x, t) + b{x)p{x, t), (32) 

where nx is the unit outer normal at a boundary point x. The probability Pn that a live 
virus arrives to a nuclear pore and the conditional mean time can be expressed using the 
SPDF as 

p poo 

Pn = Pr{r'' <t''}= (b / J{x, t) ■ nxdSx dt (33) 

JdNa Jo 
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and 

POO 

Tn = EIt" I < r'^] = / (1 - Frir" < 1 1 < r^) dt 

Jo 

POO P 

= Pn i tJ{x,t)-nxdSxdt (34) 

Jo JdNa 

where Prjr'* < t|r" < t^} is the conditional cumulative density function of the absorption 
time, given that the gene vector arrives alive to a nuclear pore. For n identical nuclear 
pores modeled as absorbing disks of radius e and when the drift is a gradient potential 
(b (x) = — V$ (aj)), the leading order term of P„ and r„ in e have been estimated |103] . 

^-^o/D 

]_ / g-$(x)/Z}^(^)^^ -$0//^ 



and, 



1 



e 



= '"^"1/" ^ — -, (36) 



L_ [ e-^(-)/Dk{x)dx + e-^o/D^ 



where $0 is the constant value of the radial potential ^{x) on the centered nucleus where 
the nuclear pores are uniformly distributed. These asymptotic results rely on the narrow 
escape theory [Slj, which is a general asymptotic method to estimate the mean first passage 
time of a Brownian motion confined in a domain to escape through a small opening. 

As the number of nuclear pores n remains small, the asymptotic validity of the previous 
expressions [361 agree in some limits with Brownian simulations |104j . However, these formulas 
do not account for the possible interactions between the small absorbing pores and indeed, 
lim„^oo,ne2<i(T) = 0. Because the nucleus contains a large number of nuclear pores (n = 2000 
[105] ) that cover a small total area (around 1% of the nuclear surface), a refined analysis has 
been developed to account for the nuclear geometry |110j . In addition, interaction between 
absorbing windows can drastically affect the MFPT |106l I107j . 

Using electrostatic arguments, when the number of nuclear pores n ^ 1 is large and 
cover uniformly a small surface of the nucleus S, the leading order of the narrow escape time 
for a pure Brownian particle is |108j 









~d' 







where Cy. is the capacity of the nucleus (for a sphere of radius 5, Cs = 47r5). When we 
neglected the geometrical interactions between the holes 1/Cs = 0, the mean time ( 137|) 
reduces to expression ( 136|) for a Brownian particle = with no killing activity k = Q. 
Interestingly, both a quantitative and qualitative information can be derived from formula 
(1371) : the MFPT is shorter for a sphere covered with many small holes compared to a single 
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large hole |104] with same surface. Actually, for a pure Brownian process, many absorbing 
small holes is equivalent to an almost totally absorbing surface jlU9j . 

Finally, in |11U] . we accounted for the interactions between absorbing nuclear pores, and 
obtained refined estimates for the probability -P„ and r„ (formula [36]) . when the velocity is 
described by the stochastic equation (!30|l 

/as \AnDe DC^JJ^ 



and 



+ 1 / ^-Hx)/Ddx 

n. ^ ^yJ" , , , . (39) 



f e-^(^)/Ddx+(-^ + -lA [ Hx)e-^(^)/Ddx 



Interestingly, for a biological cell with a spherical nucleus (radius 6 = 5fim |105] ). the n = 
2000 circular nuclear pores |105] (radius e = 25nm |105] ) cover a surface {nne^) / (47r5^) ^ 1% 
of the total nuclear surface and 1 / {4nDe) is only one third of 1/ C^,. Using the parameter for 
the Adeno-Associated- Virus, with an affective diffusion constant D = 1.3/im^s~^ [95j and 
no degradation activity {k = 0), moving in a ball of radius R = 15/im ([lllj) containing a 
nucleus and a radial potential $(r) = 0.2r/i^ms~^ |104] . using formula |39| we found that the 
time to a nuclear pore is r„ ^ Imin., which is three time longer than the estimation using 
formula fl36p . when the effect of the nucleus is not taken into account. 



6 Conclusion 

Delivering a plasmid DNA in a cell is still a challenging task and a daunting hurdle of modern 
drug delivery methods. Being able to quantify cytoplasmic trafficking cannot only be used 
to optimize gene delivery, but also to design specific drug strategy against viral infection. 

In the first part of the review, we summarized some recent progresses about quantifying 
the endosomal escape of synthetic polycations vectors, which defines a limiting step in gene 
delivery. It remains an open question to understand the precise mechanisms by which a 
plasmid is escaping and how the endosomal membrane is disrupted. We believe that any 
understanding can lead to a rational method to optimize the construction of polycations. 
In the second part, we presented some general modeling approaches to quantify cytoplasmic 
trafficking of enveloped viruses, such as influenza. It remains to find strategies to design these 
new synthetic gene vectors and for that goal, viruses appear as optimal models. Hybrid 
vectors made of a mixture of viral glycoproteins and polycations seem to be a promising 
direction and physical modeling and numerical simulations are certainly seducing tools to 
optimize the construction of such nano-scale carriers. 
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